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Abstract. We simulate anisotropic outflows of AGN, and investigate the large-scale impact 
of the cosmological population of AGN outflows over the Hubble time by performing N- 
body ACDM simulations. Using the observed quasar luminosity function to get the redshift 
and luminosity distribution, and analytical models for the outflow expansion, AGNs are 
allowed to evolve in a cosmological volume. By the present epoch, 13 - 25% of the total 
volume is found to be pervaded by AGN outflows, with 10~ 9 G magnetic field. 
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1. Introduction 

Outflows from AGN are observed in a wide va- 
riety of forms: radio galaxies, broad absorp- 
tion line quasars, Seyfert galaxies exhibiting 
intrinsic absorption in the UV, broad emis- 
sion lines, warm abs orbers and absorption 
lines in X-ra ys (e.g., ICrenshaw et al.l l2003t 
lEveretll |2007). There have been studies on 
the cosmologic al impact of quasar ou tflows 
in large scales (Furlanetto & Loeb 2001, here- 
after FL01: IScannanieco & Ohl 12004 here- 
after S O04: iLevine & Gnedinll2005[ hereafter 
LG05). iBarail (120081) investigated the cosmo- 
logical influence of radio galaxies over the 
Hubble time. All these studies considered 
spherically expanding outflows. 

On cosmological scales an outflow is ex- 
pected to move away from the high den- 
sity regions of large-scale structures, with 
the outflowing matter getting channelled 
into low-density regions of the Universe 



dMartel & Shapiroll200ll) . We implement such 
anisotropic AGN outflows within a cosmolog- 
ical volume. The simulation methodology is 
given in $2] and the results are discussed in §0 

2. The Numerical Setup 

2.1 . N-body Simulation and Distribution 

We simulate the growth of large-scale struc- 
tures in a cubic cosmological volume of 
comoving size Lbox=128/z _1 Mpc. We use 
the Particle-Mesh (PM) algorithm, with 256 3 
equal-mass particles, on a 512 3 grid. A particle 
has a mass of 1 .32 x 1O 1O M , and the grid spac- 
ing is A = 0.25/i -1 Mpc. We consider a con- 
cordance ACDM model with the cosmological 
parameters: D.m - 0.268, £2a = 0.732, Ho = 
70.410ns- 1 Mpc" 1 , Q. b = 0.0441, n s = 0.947, 
and Tcmb = 2.725. 

The redshift-dependent luminosity distri- 
bution of AGN is obtained from the bolo- 
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metric quasar lum inosity function (QLF) 
dHopkinset al.ll2007l) . 



<P(L,z) = 



4>* 



dlogL (L/L+y* +(L/L+)n 



(1) 



which gives the number of quasars per unit co- 
moving volume, per unit log 1() of luminosity. 
A fraction /outflo w = 0-2 of AGN are considere d 
to host outflows dGangulv & Brofhertorj2 008). 
The number of outflows within the simulation 
box of comoving volume Vb ox = Lj^ ox , at epoch 
z, between [L, L + dL] is, 

N{L, z) = /outflow^, zMlog 10 L] Vbox- (2) 

The AGN activity lifetime is taken as 7agn = 
10 8 yr; the maximum and minimum AGN lu- 
minosities as 10 8 L Q and 10 I4 L o . 

Using the QLF, we obtain the entire cos- 
mological population of AGN in the simulation 
volume starting from z — 6, namely the birth 
redshift (zt,ir), switch-off redshift (Zos) and 
bolometric luminosity (Lboi) of each source. A 
total of 929805 sources were produced. 

At each timestep, we filter the density 
distribution on the 512 3 grid (from the PM 
code) using a gaussian filter containing a mass 
10 10 M o , assumed as the minimum mass of a 
halo hosting an AGN. We identify the density 
peaks, or the grid cells where the filtered den- 
sity exceeds the values at the 26 neighboring 
grids. We consider the peaks that have a filtered 
density > 5 x the mean density of the box, and 
each new AGN born during that epoch (whose 
Zbir values fall within the timestep interval) is 
located at the center of one such peak cell, se- 
lected randomly. After their initial distribution, 
the AGNs are allowed to evolve according to 
the prescription in 



2.2. Outflow Model 

Despite the observational differences between 
various outflows, we stress that the AGNs host- 
ing outflows constitute a random subset of the 
whole AGN population, and we simply assume 
the same outflow model for all AGNs (also 
FL01, SO04, LG05). We allow each outflow 
to evolve through an active-AGN life ( jH2.2.1l l. 
when Zbb > z > Zos- After the central engine 



has stopped activity (when z < Zoff), it enters 
the late-expansion phase ( 32.2.21 1. 

The baryonic ambient gas density, p g (z, r), 
is considered to follow the dark matter den- 
sity, pm(z, r), in the N-body simulation: p g - 
(flblQ.u)pM- The external gas pressure is 
p g {z, r) = p g (z, r)kT g lp. The external tempera- 
ture is fixed at T g = 10 4 K assuming a photo- 
heated ambient medium, and p = 0.61 1 amu is 
the mean molecular mass. 

2.2.1. The Active Life 

The AGN activity period is short compared to 
the Hubble time. So we neglect energy losses 
and Hubble expansion of the cosmological vol- 
ume when the quasar is active. We approximate 
the shape of the outflow as spherical. 

An active outflow i s inflated by twin colli - 
mated relativistic jets (Begel man et al.l il984). 
each of length R. We consider that the ki- 
netic luminosity carried by each jet is a con- 
stant fraction of the bolometric luminosity: 
L K = eyL h „i/2, wit h e K = 0.1 (FL01; LG05; 
Shank ar et alj|20 08). The jet advance speed is 
obtained by balancing the jet momentum flux 
with the ram pressure of the ambient medium: 
Lk/(A s c) = p g (dR/dt) 2 . Here, A s is the area 
of the shocked "working" surface at the jet 
head. We use A s = 2jtR 2 8 2 , assuming that the 
shock front has a constant half-opening angle 
of 6 S = 5° relative to the central AGN (FL01). 

All the kinetic energy transported along the 
jets during an AGN's age f age = t(z) - f(Zbir) 
is transferred to the outflow. The outflow en- 
ergy is Eq = 2L K t age , and its pressure follows 
PqVq = (To - l)£b- The adiabatic index of the 
relativistic outflow plasma is To = 4/3. The 
outflow volume during this active spherical ex- 
pansion is Vq(z) = 47r/? 3 /3. 

2.2.2. The Late Anisotropic Expansion 

When AGN activity ends, the left-over high- 
pressure outflow expands into the large scales 
of the IGM with an anisotropic geometry. It 
is represented as a "bipolar spher ical cone" 
with r adius R and opening angle a dPieri et al.l 
120071) . and it follows the direction of least re- 
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sistance (DLR). We perform a second-order 
Taylor expansion of the density around each 
peak, whose coefficients are determined by 
performing a least-square fit to all the grid 
points within a distance 2A from the peak. 
We then rotate the coordinate axes such that 
the cross-terms vanish to give: p(x',y',z') = 
Ppeak - Ax' 2 - By' 2 - Cz' 2 - The largest of the 
coefficients A, B, C gives the DLR. 

The Sedov-Taylor adiabatic blast wave 
model is used to ob tain the radius of th e over- 
pressured outflow (ICastor et al.l Il975b SO04; 

LG05): R(z) = & (E to ttl g Jp- g (z)) U5 . We obtain 
p^(z) by averaging the gas density of the grid 
cells occurring within the outflow volume. For 
a strong explosion in the F g = 5/3 ambient 
gas, = 112. Here the total kinetic energy 
injected into the outflow by the AGN is E tot = 
2LaTagn. Adiabatic expansion losses are con- 
sidered, and the outflow pressure evolves as 
poR 3rs = constant, with the constant derived 
from the values at the end of the active phase. 

The outflow follows an anisotropic expan- 
sion as long as its pressure exceeds the exter- 
nal pressure of the IGM, i.e., po(z) > p g (z). 
During this late biconical expansion, Vq(z) = 
4nR\l - cos(or/2))/3. 

When po{z) < p g {z), or, the outflow has 
reached pressure equilibrium with the IGM, it 
has no further expansion. After this, the out- 
flow simply evolves passively with the Hubble 
flow of the cosmological volume. Thus an out- 
flow in pressure equilibrium attains a final vol- 
ume of V = 4ttR 3 JI - cos(o'/2))/3, where R f 
is the final comoving radius of the outflow. 



3. Results and Discussion 

At each timestep the total volume occu- 
pied by the AGN outflows is computed by 
counting the contributions of all the sources 
born by then, both the active ones and those 
in the anisotropic phase. We performed 4 
simulations with opening angles of a = 
60°, 90°, 120°, 180°, all with e K = 0.1, and one 
with a = 120° and €k - 0.05, whose results 
are shown in the figures. 

We count the grid cells in the simulation 
box which occur inside the volume of one or 
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Fig. 1. Volume of the simulation box filled by AGN 
outflows (Wagn) a s a fraction of the total volume 
(solid), and as a fraction of volumes of various over- 
densities: N(p > p) (dash dot), N(p > 2p) (dashed), 
N(p > 3p) (dotted), N(p > Sp) (dash dot dot dot), 
N(p > 7p) (long dashes). 



more AGN outflows. The total number of these 
filled cells, Nagn, gives the total volume of the 
box occupied by outflows. We express the total 
volume filled as a fraction of volumes of vari- 
ous overdensities in the box, N p = Nip > Cp), 
where p = (1 + z) 3 Qm3//q/(87tG) is the mean 
matter density of a spatially flat Universe (the 
box) at an epoch z- So N p gives the number of 
cells which are at a density C times the mean 
density. We find N p for C = 0,1,2,3,5,7; 
C — gives the total volume of the box. 

Figure[T|shows the redshift evolution of the 
volume filling factors for our 5 simulation runs. 
With 10% kinetic efficiency, 0.13 of the entire 
Universe is filled at present by AGN outflows 
with an opening angle of 60°; the fraction in- 
creases to 0.17 with 90°, 0.21 with 120°, and 
0.25 with 180°. A 5% kinetic efficiency and 
a = 120° fills 0.13 of the volume. In all our 
runs, the outflows fill up all of the regions with 
p > 2p by z = 0.3. With e K = 0.1 and a = 90° 
or higher, the outflows permeate all the over- 
dense regions (p > p) by z = 0.1. 
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ates the filled ov erdense volumes , consistent 
with the results of lRyu et ai1 (l2008). At a given 
redshift, the energy density and magnetic field 
are larger for smaller opening angles of the 
anisotropic outflows. 

We conclude that, using our N-body sim- 
ulations, the cosmological population of AGN 
outflows pervade 13 - 25% of the volume of 
the Universe by the present. A magnetic field 
of ~ 10~ 9 G is infused in the filled volumes at 
z = 0. 

Acknowledgements. All calculations were per- 
formed at the Laboratoire d'astrophysique 
numerique, Universite Laval. We thank the Canada 
Research Chair program and NSERC for support. 



Fig. 2. Volume fraction filled by AGN outflows 
(top), the volume weighted average of the total en- 
ergy density inside outflow volumes (u E ) (middle), 
and the equipartition magnetic field within the filled 
volumes (So) (bottom). 



It is the overdense regions of the Universe 
which gravitationally collapse to form stars 
and galaxies. So evidently the AGN outflows 
have a profound cosmological impact on the 
protogalactic regions. We note that the vol- 
umes obtained by LG05 (100% filling by z ~ 
1) are higher than ours. 

We perform preliminary estimates of the 
energy density and magnetic field in the vol- 
umes of the Universe filled by the AGN out- 
flows. The energy density inside the outflow 
behaves similar to the outflow pressure evolv- 
ing adiabatically ( jj2.2.21 >. ue = 3po. Assuming 
equipartition of energy between magnetic field 
of strength Bq and relativistic particles inside 
the outflow, the magnetic energy density is 
ub - ue/2 = Bq/(8jt). We define the volume 
weighted average of a physical quantity J?l as 
(&)(z) = 2G?lVu)/ 2 V , where the summa- 
tion is over all outflows existing in the simula- 
tion box at that epoch. 

Figure [2] shows the redshift evolution of 
the total volume filling fraction, (ue) and (Bo). 
The energy density and magnetic field decrease 
with redshift as larger volumes are filled. At 
z — 0, a magnetic field of ~ 10~ 9 G perme- 
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